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Abstract. We present high resolution echelle spectra taken during four observing runs from 1999 to 2001 of 
the recently X-ray/EUV selected chromospherically active binary BK Psc (2RE J0039-I-103). Our observations 
confirm the single-lined spectroscopic binary (SBl) nature of this system and allow us to obtain, for the first time, 
the orbital solution of the system as in the case of a SB2 system. We have determined precise radial velocities 
of both components: for the primary by using the cross correlation technique, and for the secondary by using its 
chromospheric emission lines. We have obtained a circular orbit with an orbital period of 2.1663 days, very close 
to its photometric period of 2.24 days (indicating synchronous rotation). The spectral type (K5V) we determined 
for our spectra and the mass ratio (1.8) and minimum masses (Msin'^i) resulting from the orbital solution are 
compatible with the observed K5V primary and an unseen M3V secondary. Using this spectral classification, the 
projected rotational velocity (wsini, of 17.1 km s~^) obtained from the width of the cross-correlation function 
and the data provided by Hipparcos, we have derived other fundamental stellar parameters. The kinematics 
and the non- detection of the Li I line indicate that it is an old star. The analysis of the optical chromospheric 
activity indicators from the Ca ll H & K to Ca ll IRT lines, by using the spectral subtraction technique, indicates 
that both components of the binary system show high levels of chromospheric activity. Ha emission above the 
continuum from both components is a persistent feature of this system during the period 1999 to 2001 of our 
observations as well as in previous observations. The Ha and H/3 emission seems to arise from prominence-like 
material, and the Ca ll IRT emission from plage-like regions. 

Key words, stars: individual: BK Psc - stars: activity - stars: binaries: spectroscopic - stars: chromospheres - 
stars: late-type - stars: rotation 



1. Introduction 

This paper is a part of our ongoing project on multiwave- 
length optical observations aimed at studying the chromo- 
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sphere of active binary systems. For this purpose we use 
the information provided by several optical spectroscopic 
features that are formed at different heights in the chro- 
mosphere (see Montes et al. 1997, Paper I; Montes et al. 
1998, Paper II; Montes et aL 2000, Paper III). In addition 
to study stellar activity, the high resolution spectroscopic 
observations we use in this project allow us to determine 
radial velocities and to obtain and improve fundamental 
stellar parameters. 

In this paper we focus our attention on the X-ray/EUV 
selected chromospherically active binary BK Psc (2RE 
J003939+103925, BD-l-09 73, G 1-10, LHS 1118). It is 
a high proper-motion star with photometry reported by 
Stephenson (1986), Sandage & Kowal (1986) and Weis 
(1991) (V^=10.5; U -B = 0.97; B -V ^ 1.17; V - R = 
0.73; R-I = 0.60). Bidelman (1985) gives a K5 spectral 
type for this star that was confirmed later by Jeffries et 
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Table 1. Observing log 
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al. (1995), but Stephenson (1986) listed it as K4:p. Due 
to a blue excess in the {U — B) color index the presence 
of a white dwarf companion has been suggested by Weis 
(1991) and Cutispoto et al. (1999). Radial velocity vari- 
ations indicate that it is a binary system (Jeffries et al. 
1995; Cutispoto et al. 1999), but no orbital solution has 
been determined until now. BK Psc was detected as an 
extreme ultraviolet (EUV) source by the ROSAT Wide 
Field Camera all-sky survey (Pounds et al. 1993; Pye et 
al. 1995). The ehromospheric activity of this star was de- 
tected in the optical identification program of ROSAT 
EUV Sources by Mason et al. (1995) {EW {Y{a)=l.l A and 
EW{C&ii K)=2.7 A) and Jeffries et al. (1995) (strong Ha 
emission above the continuum, £'T^(Ha)=1.0 A). Finally, 
Cutispoto et al. (1999) found that the best fit for their ob- 
served colors (y=10.43; U-B = 0.92; B-V = 1.16; V-R 
= 0.73; V -I = 1.49) and HiPPARCOS distance (d = 32.8 
pc) is obtained by assuming a K5V or K6V primary, a 
M4V secondary and a possible white dwarf (WD). These 
authors also confirm the optical variability of this star 
(photometric period Pphot = 2.24 ±0.04 days), previously 
reported in the SAAO Annual report (1993) (Pphot = 2.17 
days). 

In this paper we present high resolution echelle spectra 
of BK Psc, obtained at different epochs, that allow us to 
measure radial and rotational velocities by using the cross- 
correlation technique. With these observations we confirm 
the binary nature (SBl) of this system and determine its 
orbital solution for the first time. We have obtained an 
orbital period of 2.17 days, very close to the photometric 
period of 2.24 days, indicating nearly synchronous rota- 
tion. Furthermore, we have applied the spectral subtrac- 
tion technique to study the ehromospheric excess emission 
in the Ca ii H & K, Ca ii IRT, Ha and other Balmer lines 
from the primary and secondary components of the sys- 
tem. Preliminary results for this system can be found in 
Galvez et al. (2001) and Montes et al. (2001a). 

In Sect. 2 we give the details of our observations and 
data reduction. In Sect. 3 the procedures to obtain the 
stellar parameters and the orbital determination of the 
binary system are described in more detail and the re- 
sults are discussed. The individual behavior of the dif- 
ferent ehromospheric activity indicators is described in 
Sect. 4. Finally in Sect. 5 we give the conclusions. 



2. Observations and data reduction 

The spectroscopic observations of BK Psc were obtained 
during four observing runs: 

1) 2.2 m-FOCES 1999/07 

This run took place on 24-29 July 1999 using the 
2.2 m telescope at the German Spanish Astronomical 
Observatory (CAHA) (Almeria, Spain). The Fibre Optics 
Cassegrain Echelle Spectrograph (FOCES) (Pfeiffer et al. 
1998) was used with a 2048x2048 LORAL#lli CCD 
detector. The wavelength range covers from 3910 to 9075 
A in 84 orders. The reciprocal dispersion ranges from 0.03 
to 0.07 A/pixel and the spectral resolution, determined 
as the full width at half maximum (FWHM) of the arc 
comparison lines, ranges from 0.09 to 0.26 A. 

2) INT-MUSICOS 2000/08 

This run was done on 11-14 August 2000, with the 
2.5 m Isaac Newton telescope (INT) at the Observatorio 
del Roque de Los Muchachos (La Palma, Spain) using 
the ESA MUSICOS Spectrograph. This is a fibre-fed 
cross-dispersed echelle spectrograph (Baudrand & Bohm 
1992) and developed as a part of MUlti-SIte Continuous 
Spectroscopy (MUSICOS) project. During this observing 
run, a 2148x4148 pixel EEVlOa (EEV4280) CCD detector 
was used. The wavelength range covers from 4000 to 10000 
A in 87 orders. The reciprocal dispersion ranges from 0.03 
to 0.08 A/pixel and the spectral resolution (FWHM) from 
0.06 to 0.24 A. 

3) NOT-SOFIN 2000/11 

This run took place on 6-9 November 2000 using the 
2.56 m Nordic Optical Telescope (NOT). The Soviet 
Finnish High Resolution Echelle Spectrograph (SOFIN) 
was used with an echelle grating (79 grooves/mm), cam- 
era ASTROMED-3200 and a 1152x298 pixel EEV P88100 
CCD detector. The wavelength range covers from 3900 to 
8400 A in 37 orders. The reciprocal dispersion ranges from 
0.06 to 0.14 A/pixel and the spectral resolution (FWHM) 
from 0.15 to 0.42 A. 

4) 2.2 m-FOCES 2001/09 

The last run took place on 21-24 September 2001 using 
the 2.2 m telescope at the German Spanish Astronomical 
Observatory (CAHA) (Almerfa, Spain). In this run, the 
FOCES spectrograph (Pfeiffer et al. 1998) was used with 
a 2048x2048 (24/i) SITE#ld CCD detector. The wave- 
length range covers from 3450 to 10700 A in 112 orders. 
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Table 2. Stellar parameters of BK Psc 
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Tgp 


SB 


V 


B-V V-R Pi,b 
(days) 


fphot 

(days) 


V sin 
(km s"^) 


TT 

(mas) 


/ia cos 5 

(mas/yr) 


(mas/yr) 


K5/6:V/M4:V + WD 


1 


10.43 


1.16 0.73 2.1663 


2.24 


17.1 


30.52±1.79 


524.9±2.1 


-198.0±2.1 



^ values determined in this paper 



The reciprocal dispersion ranges from 0.04 to 0.13 A/pixel 
and the spectral resolution (FWHM) ranges from 0.1 to 
0.34 A. 

In Table |l| we give the observing log. For each obser- 
vation we list date (Day), universal time (UT), exposure 
time (Exp), and the signal to noise ratio (S/N) measured 
in the Ha line region. 

The spectra have been extracted using the standard 
reduction procedures in the IRAF|^ package. The wave- 
length calibration was obtained by taking spectra of a 
Th-Ar lamp. Finally, the spectra have been normalized 
by a polynomial fit to the observed continuum. 

3. Stellar parameters of the binary system 

The adopted stellar parameters of BK Psc are given in 
Table ^. Spectral type and the photometric data {V , 
B — V, V — R, Pphot) are taken from Cutispoto et al. 
(1999). Orbital period (Porb) and projected rotational ve- 
locity (wsini) have been determined in this paper (see 
below). The astrometric data (parallax, it, proper mo- 
tions, fiaCOsS and ^s) are from Hipparcos (ESA 1997) 
and Tycho-2 (H0g et al. 2000) catalogues. 

3.1. Spectral types 

In order to obtain the spectral type of this binary sytem 
we have compared our high resolution echelle spectra, in 
several spectral orders free of lines sensitive to chromo- 
spheric activity, with spectra of inactive reference stars 
of different spectral types and luminosity classes observed 
during the same observing run. This analysis made use of 
the program STARMOD developed at Penn State University 
(Harden 1985) and modified later by us. With this pro- 
gram a synthesized stellar spectrum is constructed from 
artificially rotationally broadened, radial-velocity shifted, 
and weighted spectra of appropiate reference stars. We 
have obtained the best fit between observed and synthetic 
spectra when we use a K5V primary component without 
any contribution from a secondary component. As refer- 
ence stars we have used the K5V stars HD 154363 for the 
first run, 61 Cyg A for for the other three runs. This spec- 
tral classification is in agreement with K5 spectral type 

^ TRAP is distributed by the National Optical Observatory, 
which is operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation. 



reported by Hidelman (1985) and Jeffries et al. (1995), 
and the K5/6:V + M4:V classification given by Cutispoto 
et al. (1999), where the M4V secondary has no apreciable 
contribution to the spectra. 

3.2. Rotational velocity 

The projected rotational velocity (i>sini) of this star has 
been previously estimated as < 12 km s~^ (Jeffries et al. 
1995) and 18±2 km s^^ by Cutisposto et al. (1999). 

Hy using the program STARMOD we have obtained the 
best fits for each observing run, with v sin i values around 
15-18 km s^^. In order to determine a more accurate ro- 
tational velocity of HK Psc we have made use of the cross- 
correlation technique in our high resolution echelle spectra 
by using the routine FXCOR in IRAF. When a stellar spec- 
trum with rotationally broadened lines is cross-correlated 
against a narrow-lined spectrum, the width of the cross- 
correlation function (CCF) is sensitive to the amount of 
rotational broadening of the first spectrum. Thus, by mea- 
suring this width, one can obtain a measurement of the 
rotational velocity of the star. 

The observed spectra of HK Psc were cross-correlated 
against the spectrum of a template star (a slowly rotating 
star of similar spectral type) and the width (FWHM) of 
CCF determined. The calibration of this width to yield 
an estimation of v sin i is determined by cross-correlating 
artificially broadened spectra of the template star with 
the original template star spectrum. The broadened spec- 
tra were created for v sin i spanning the expected range of 
values by convolution with a theoretical rotational profile 
(Gray 1992) using the program STARMOD. The resultant 
relationship between vsini and FWHM of the CCF was 
fitted with a fourth-order polynomial. We have tested this 
method with stars of known rotational velocity, obtaining 
a good agreement. The uncertainties on the vsini values 
obtained by this method have been calculated using the 
parameter R defined by Tonry & Davis (1979) as the ratio 
of the CCF height to the rms antisymmetric component. 
This parameter is computed by the task FXCOR and pro- 
vides a measure of the signal to noise ratio of the CCF. 
Tonry & Davis (1979) showed that errors in the FWHM 
of the CCF are proportional to (1 -I- i?)~^ and Hartmann 
et al. (1986) and Rhode et al. (2001) found that the quan- 
tity ±v s'mi{l + R)~^ provide a good estimate for the 90% 
confidence level of a vsini measurement. Thus, we have 
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Orbit of BK Psc 
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Fig. 1. Radial velocity data and fit vs the orbital phase. Solid circles represent the primary and open circles represent 
the secondary. The solid curves represent a minimum fit orbit solution as described in the text 



adopted ±usini(l + R) ^ as a reasonable estimate of the 
uncertainties on our wsini measurements. 

As template stars for BK Psc we have used the slowly 
rotating K5V stars above mentioned in the spectral type 
classification. We have determined v sin i by this method 
in all the spectra of BK Psc available and the result- 
ing weighted means in each observing run are 16. 2 ±2.1, 
16.4±0.7, 19.5±0.8, and 15.0±1.0km s'^ respectively. The 
weighted mean for all the observing runs is 17.1±0.5 km 
s""'^, which is the value given in Table 0. 



3.3. Radial velocities and orbital parameters 

Recently, Jeffries ct al. (1995) and Cutispoto ct al. (1999) 
have noted that BK Psc is a single-lined spectroscopic bi- 
nary (SBl), but they provided only four radial velocity 
measurements. The detailed analysis of our spectra and 
our radial velocities measurements, confirm the SBl na- 
ture of this system. In our spectra only the photospheric 
absorption lines coming from the primary component are 
observed throughout all the spectral range. On the con- 
trary, the chromospheric emission lines from both compo- 
nents are detected in our spectra (see Fig. |2|), and thereby 
it has been possible to measure the radial velocity of the 
secondary and in this way to obtain the orbital solution of 
the system as in the case of a double-lined spectroscopic 
binary (SB2). 



For the primary component of BK Psc the heliocentric 
radial velocities have been determined by using the cross- 
correlation technique. The spectra of BK Psc were cross- 
correlated order by order, by using the routine FXCOR in 
IRAF, against spectra of radial velocity standards with 
similar spectral type taken from Beavers et al. (1979). The 
radial velocity was derived for each order from the position 
of the cross-correlation peak, and the uncertainties were 
calculated by fxcor based on the fitted peak height and 
the antisymmetric noise as described by Tonry & Davis 
(1979). In Table ^ we list, for each spectrum, the helio- 
centric radial velocities (T4iei) and their associated errors 
((Ty ) obtained as weighted means of the individual values 
deduced for each order. Those orders that contain chro- 
mospheric features and prominent telluric lines have been 
excluded when determining the mean velocity. 

For the secondary component, we have used the infor- 
mation provided by the chromospheric emissions that are 
detected for both components in the Ha, Ca ii H & K and 
other Balmer lines. The contribution of each component 
to the observed profile has been deblended by mean of a 
two Gaussian fit (see Fig. ||) and the relative wavelength 
separation of the secondary component with respect to the 
primary has been used to determine its heliocentric radial 
velocity (listed in Table ||) . 

We have computed the orbital solution of BK Psc us- 
ing our eleven data of radial velocities (for both compo- 
nent) and the four values given (only for the primary) 
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Table 3. Radial velocities of BK Psc Table 4. Orbital solution of BK Psc 
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Obs. 



HJD 

2400000+ 



Primary 



Vhcl 

(km s 



± cry 
1\ 



Secondary 

Vhcl ± (TV 

(km s"^) 



1992 (J95)^ 
1992 (J95)^ 

1995 (C99)^ 
1995 {C99f 



48845.581 
48851.540 

49956.718 
49958.774 



2.9 ± 2.0 
-56.9 ± 2.0 

-26.3 ± 4.0 
-40.6 ± 4.0 



2.2 m 1999 


51385, 


,615 


-29, 


,80 


± 


0, 


,29 


31.92±4.3 


2.2 m 1999 


51387, 


,650 


-9, 


.67 


± 


0, 


,46 


-9.67±4.5 


2.2 m 1999 


51389, 


,633 


16, 


,51 


± 


0, 


,47 


-57.65±4.5 


INT 2000 


51767, 


,692 


-29, 


,98 


± 


0, 


,13 


36.79±4.1 


INT 2000 


51768, 


,689 


22, 


,71 


± 


0, 


,11 


-70.23±4.1 


INT 2000 


51769, 


,706 


-52, 


,42 


± 


0, 


,13 


55.93±4.4 


INT 2000 


51770, 


,703 


37, 


,91 


± 


0, 


,09 


-93.27±4.1 


NOT 2000 


51854, 


,551 


-6, 


.84 


± 


0, 


,21 


-6.84±4.2 


NOT 2000 


51855, 


,470 


3, 


,06 


± 


0, 


,17 


-46.55±4.2 


NOT 2000 


51856, 


,527 


-32, 


,54 


± 


0, 


,14 


31.00±4.1 


NOT 2000 


51857, 


,486 


27, 


,93 


± 


0, 


,78 


-75.80±4.8 


2.2 m 2001 


52176, 


,554 


-56, 


,65 


± 


0, 


,23 


71.82±4.2 


2.2 m 2001 


52177, 


,448 


16, 


,25 


± 


0, 


,10 


-57.80±4.1 



1 J95: Jeffries et al. (1995), 
^ C99: Cutispoto et al. (1999). 



by Jeffries et al. (1995) and Cutispoto et al. (1999) (see 
Tabled) . The radial velocity data (Table |^) are plotted in 
Fig. ^Solid circles represent the primary and open cir- 
cles represent the secondary. The curves represent a min- 
imum fit orbit solution. The orbit fitting code uses 
the Numerical Recipes implementation of the Levenberg- 
Marquardt method of fitting a non-linear function to the 
data and weights each datum according to its associated 
uncertainty (Press et al. 1986). The program simultane- 
ously solves for the orbital period, Forb, the epoch of pe- 
riastron passage, Tq, the longitude of pcriastron, uj, the 
eccentricity, e, the primary star's radial velocity ampli- 
tude, Ki, the heliocentric center of mass velocity, 7, and 
the mass ratio, q. The secondary star's radial velocity am- 
plitude, K2, is qKi. The orbital solution and relevant de- 
rived quantities are given in Table |. In this Table we give 
Tconj as the heliocentric Julian date on conjunction with 
the hotter star behind, in order to adopt the same crite- 
ria used by Strassmeier et al. (1993) in their catalog of 
chromospherically active binary stars. We have used this 
criterion to calculate the orbital phases of all the observa- 
tions reported in this paper. 

This binary system results in a circular orbit (e = 
0.0025) with an orbital period of 2.1663 days, which is 
very similar to its rotational period derived from the pho- 
tometry (Pphot = 2.24 days) indicating nearly synchronous 
rotation. 



Element 


Value 


Uncertainty 


Units 


P K 

orb 


2.1663 


0.0015 


days 


^conj 


2451383.32 


0.20 


HJD 


UJ 


84.88 


0.61 


degrees 


g 


0025 


0.0074 






52.70 


0.73 


km s ^ 


K-, 


95.09 


2.66 


km s~ ^ 


'V 

/ 


-10.95 


0.32 


km s ^ 


q = M1/M2 


1.80 


0.04 






1.57 


0.02 


10® km 




2.83 


0.08 


10^ km 


d sixii 


4.40 


0.08 


lO'' km 




0.0294 




AU 


5? 


6.33 




R.0 


Ml sin^i 


0.466 


0.027 




M2 sin^i 


0.258 


0.015 


Me 


/(M)i 


0.0328 


0.0013 


Me 


/(M)2 


0.1917 


0.0025 


Me 



3.4. Other derived quantities 

For the K5V spectral type of the observed primary 
component we can adopt from Landolt-Bornstein tables 
(Schmidt-Kaler 1982) a primary mass Mi — 0.67 Mq. 
According to the mass ratio {q — M1/M2 = 1.80) from 
the orbital solution we estimate for the secondary a mass 
M2 = 0.37 Mq which corresponds (Schmidt-Kaler 1982) 
to a M3V star. With these spectral types the difference in 
bolomctric magnitudes between both components is 2.5 
and the difference in visual magnitudes is 3.9 which is in 
agreement with an unseen secondary component and the 
spectral classification reported by Cutispoto et al. (1999). 

We have estimated the radius of the primary compo- 
nent by using the parallax (30.52 mas) given by Hipparcos 
(ESA 1997) and the unspotted V magnitude, taken as 
the brightest magnitude (10.43) of the values given by 
Cutispoto et al. (1999). This V magnitude is very close 
to the value given by Hipparcos {Vr — 10.60 that corre- 
sponds to V — 10.48). As the system is relatively close, 
to calculate the absolute magnitude My, no interstel- 
lar reddening was assumed. The bolomctric correction 
{BC ~ —0.72) corresponding to the K5V primary from 
(Schmidt-Kaler 1982) has been used to compute the bolo- 
mctric magnitude, Mboi and luminosity, L/Lq. Assuming 
that the contribution of the secondary to this total lu- 
minosity is very small, we have used this L/Lq and the 
effective temperature {Tf.ff = 4350 K) corresponding to 
a K5V to determine a radius (we called Rhip) for the 
primary Ruip — 0.60 ± 0.04i?Q. The errors in these de- 
rived quantities are dominated by the error in the par- 
allax (±1.79 mas) given by Hipparcos and in the Te// 
(±100 K). This radius can be compared with an inde- 
pendent determination of the minimum radius (i?sini). 
Taking as rotational period the photometric period (2.24 
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days) given by Cutispoto et al. (1999) and the rota- 
tional velocity vsini — 17.1 determined by us we found 
Rsini = 0.76±0.03i?Q. The error in this case is dominated 
by the uncertainty in vsini. This value of Rsini should 
be smaller than Rnip, but we have found a value slightly 
larger. Within the errors, however, the agreement between 
both radii is acceptable. In addition, the radius for a K5V 
in Schmidt-Kaler (1982), R = 0.727?©, is halfway between 
Rnip and Rsini. This low value of Rnip also suggests 
that the V magnitude we have used can be effected by cool 
spots on the stellar surface. Using the measured minimum 
radius i?sini and the effective temperature (Teff = 4350 
K) corresponding to a K5V we obtain a low limit of the 
of the stellar luminosity {L — 0.167Lq) and brightness 
{V — 9.88). Using this luminosity and the the mass- 
luminosity relation for main sequence stars we obtain an 
estimate of the mass of the primary (0.594 Mq) that is 
compatible with the mass of a K5V-K7V. With this mass 
and the minimum mass obtained from the orbital solution 
the mass for the secondary is M2 = 0.33 Mq which corre- 
spond to a M3V star, similar to the result obtained with 
the first method described above. 

We can summarize the adopted and derived quantities 
of the primary and secondary components of BK Psc as 
follows: 

f M^K5V) = 0.67, Afi™„ « 0.59 Mq 
Primary K5V = < Ri{K5V) = 0.72, Ri sini « 0.76 Rq 
[ T^K5V) - 4350 K 

( M2 « 0.37, M2„„„ w 0.33 Mq 
Secondary M3U = < R2{M3V) = 0-45 Rq 
[ T^iMzv) = 3470 K 

The inclination of the system isi — 62.4° if we compare 
the minimum mass of the primary {Mi sin'^z=0.466 Mq) 
deduced from the orbit with the mass adopted for a K5V. 
The minimum inclination angle for eclipses to occur (imm) 
is given by cosimin — (i?i+i?2)/a. Using the radii adopted 
for the primary and secondary components and the semi- 
mayor axis of the orbit (a sini) lead to imin = 80.1°. Since 
the photometric observations show no evidence of eclipses, 
the inclination of BK Psc must be lower than 80.1°, which 
is in agreement with the previous estimation of i. 

3.5. Kinematics and age 

BK Psc is a high proper-motion star included in the stud- 
ies of Stephenson (1986), Sandage & Kowal (1986) and 
Weis (1991). It is a relatively nearby star (d = 32.8 pc) 
with astrometric data measured by Hipparcos (ESA 1997) 
and Tycho-2 (H0g et al. 2000) catalogues (see Table |). 

We have computed the galactic space- velocity compo- 
nents ([/, V, W) using as radial velocity the center of mass 
velocity (7) (for details see Monies et al. 2001b). The re- 
sulting values and associated errors are given in Table |. 

The large total velocity Upotai and the U, V, W ve- 
locity components that lie clearly outside the young disk 
population boundaries in the {U, V) and {U, W) diagrams 



Table 5. Galactic space- velocity components 



U±au 


V±av 


W ± aw 


Vrotal 


(km s~^) 


(km s-i) 


(km s~i) 


(km s-i) 


-55.23 ± 3.38 


-63.27 ± 3.35 


-25.66 ± 2.04 


87.82 



(Eggen 1984, 1989; Monies et al. 2001b) indicate that BK 
Psc is an old disk star. 

The spectral region of the resonance doublet of Li i at 
A6708 A is included in all our spectra of BK Psc. The de- 
tailed analysis of the spectra indicates that this line is not 
detected in this star. As it is well known, this spectroscopic 
feature is an important diagnostic of age in late- type stars. 
In addition, it is also known that a large number of chro- 
mospherically active binaries shows Li i abundances higher 
than other stars of the same mass and evolutionary stage 
(see Paper II, III and references therein). Therefore, non- 
detection of the Li I line in this active star indicates that 
it is an old star which is in agreement with its kinematics. 

4. Chromospheric activity indicators 

The echelle spectra analysed in this work allow us to study 
the behaviour of the different indicators from the Ca 11 H 
& K to the Ca 11 IRT lines, which are formed at different 
atmospheric heights. The chromospheric contribution in 
these features has been determined by using the spectral 
subtraction technique described in detail by Montes et al. 
(1995a) and Paper I, II, and III. The synthesized spectrum 
was constructed using the program STARMOD. Taken into 
account the stellar parameters derived in Sect. 3 we have 
used only a K5V primary component without any contri- 
bution from a secondary component. The inactive K5V 
stars used as reference stars arc HD 154363 for the first 
observing run, and 61 Cyg A for the other three runs. 

In Table ^ we give the excess emission equivalent width 
(EW) (measured in the subtracted spectra) for the Ca 11 
H & K, He, m, H7, H/3, Ha, and Ca 11 IRT (A8498, A8542, 
A8662) lines. When the emission features from both com- 
ponents can be deblended, we give the EW for the hot 
and cool (H/C) components. The uncertainties in the mea- 
sured EW were estimated taking into account: a) the typ- 
ical internal precisions of STARMOD (0.5 - 2 km s^^ in ve- 
locity shifts, and ±5 km s^^ in v sini), b) the rms obtained 
in the fit between observed and synthesized spectra in the 
spectral regions outside the chromospheric features (typi- 
cally in the range 0.01-0.03) and c) the standard deviations 
resulting in the EW measurements. The final estimated 
errors are in the range 10-20%. 

The measured EWs given in Table ^ have been cor- 
rected for the relative contribution of each component to 
the total continuum determined by means of the radii and 
temperatures assumed in Sect. 3. For instance, in the Ha 
line region the relative contributions arc Sh ~ 0.94 for the 
hot component and Sc = 0.06 for the cool component, and 
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Fig. 2. Spectra of BK Psc in the Ha line region. The observed spectrum (solid-line) and the synthesized spectrum 
(dashed-line) are plotted in the left panel and the subtracted spectrum (dotted line) in the right panel. The position 
of the Ha line for the primary (P) and secondary (S) components are marked with short vertical lines. We have 
superposed the two-Gaussian fit used to deblend, in the subtracted spectrum, the excess emission coming from both 
components. 



8 



M.C. Galvez et al.: The X-ray/EUV selected binary BK Psc 



the corrected EWs for the hot and cool components are 
obtained multiplying by a factor 1/5h and 1/Sc, respec- 
tively. Finally, these corrected EWs have been converted 
to absolute surface fluxes by using the empirical stellar 
flux scales calibrated by Hall (1996) as a function of the 
star colour index. In our case, we have used the B — V 
index and the corresponding coefBcients for Ca ii H & K, 
Ha and Ca ii IRT, using the same as Ca ii H & K for 
He, and derived the H5, H7 and H/3 fluxes by making an 
interpolation between the values of Ca 11 H & K and Ha. 
The logarithm of the obtained absolute flux at the stel- 
lar surface (logFs ) for the different chromospheric activity 
indicators is given in Table |^. 

In Figs. 1^ and H we have plotted for each observation 
in the Ha and Ca 11 IRT A8498, A8542 line region the 
observed spectrum (solid-line) and the synthesized spec- 
trum (dashed-line) in the left panel, and the subtracted 
spectrum (dotted line), in the right panel. The observing 
run and the orbital phase {(p) of each spectrum are also 
given in these figures. The observed spectra in the Ca 11 
H & K line region are plotted in Fig.^, and representative 
subtracted spectra in the H/3, H7 and US line regions are 
plotted in Fig. ||. 

4.1. The Ha line 

The Ha line region is included in our spectra in the four 
observing runs. In all cases we have detected, in the ob- 
served spectra (see Fig. ^left panel), strong Ha emission 
above the continuum coming from the primary component 
and a small Ha emission coming from the secondary com- 
ponent. In all the spectra, except two which are very close 
to conjunction, we were able to deblend the emission com- 
ing from both components by using a two-Gaussian fit to 
the subtracted spectra (see Fig. ^ right panel) . 

The Ha emission of the primary exhibits a central self- 
absorption similar to that observed in many M active stars 
(Stauffer & Hartmann 1986) and some K active stars like 
the dK5e binary V833 Tau (Montes et. al. 1995b) and the 
K4V single star V834 Tau (Montes et. al. 2001c). This 
self-absorption feature is a consequence of the line for- 
mation process in the chromosphere of very active stars 
(Houdebine & Doyle 1994). 

The Ha emission of BK Psc is persistent during the 
period of time covered by our observations (from 1999 
to 2001). In addition, strong Ha emission above the con- 
tinuum from the primary component was also detected 
in previous spectra of this system taken in 1992 with 
EW{Ra) = 1.0 A (Jeffries et al. 1995) and EW(Ra)^l.l 
A (Mason et al. 1995). These EWs are lower than the 
EWs determined by us because these authors determined 
the EWs in the observed spectra and our EWs have been 
measured in the subtracted spectra, after eliminating the 
photospheric contribution. This persistent Ha emission 
detected in BK Psc indicates that it is a very chromo- 
spherically active binary system similar to some RS CVn 
systems like V711 Tau, UX Ari, HU Vir, and DM UMa, 



H-y 





Wavelength (A) Wavelength (A) Wavelength (A) 

Fig. 3. Subtracted spectra in the region of the H/3, H7, 
and H^ lines. Clear excess emission from the primary (P) 
and secondary (S) components is detected 



and some BY Dra systems like BY Dra itself, and YY 
Gem, which always show Ha emission above the contin- 
uum. 

The detection of Ha emission from the cool secondary 
component (M3V) of BK Psc indicates that this star has 
a very strong Ha emission, since its photospheric contri- 
bution to the observed continuum is practically negligi- 
ble. Strong Ha emission is typical of the group of M type 
stars called dMe stars, some of which also show a scaled- 
up version of solar flares and are known as flare stars of 
UV Get type stars. These latter stars are characterized 
by dramatic increases in the Hydrogen Balmer emission 
lines. However, this seems not to be the case for the sec- 
ondary component of BK Psc since the Ha emission we 
have detected has a similar intensity in the four observing 
runs. 

The Ha emission of the primary component shows 
small variations with the orbital phase, for instance, in 
the first run the EW changes from 1.7 to 1.1 A. Seasonal 
variations are also detected, with larger EW(Jia) in 1999 
than in 2001. 



4.2. The Hp, Hj and H6 lines 

The other three Balmer lines included in all our spectra 
(H/3, H7 and US) also show evidence of chromospheric ac- 
tivity. After applying the spectral subtraction, clear excess 
emission from both components is detected (see Fig. |[). 
When the S/N is high enough we have deblended the 
emission coming from both components by using a two- 
Gaussian fit to the subtracted spectra (see Table |^). The 
three lines show small seasonal and orbital phase varia- 
tions with the same trend that the Ha line. 

We have also measured the ratio of excess emission 
EW in the Ha and H/3 lines, , and the ratio of 

excess emission with the correction: 

^^^^i^-;. 0.2444. 2.512(^-^) 
EW{llf3) 
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given by Hall & Ramsey (1992) that takes into account 
the absolute flux density in these lines and the color dif- 
ference in the components. We have obtained for the pri- 
mary component in the range of 3 to 4 in all our 
spectra. These values indicate, according to Buzasi (1989) 
and Hall & Ramsey (1992), the presence of prominence- 
like material at the stellar surface. 

4.3. The Caw H & K and He lines 

The Ca ii H & K line region is included in the spectra 
of the FOCES 1999 and 2001 observing runs. Only the 
Ca II K line is included in the NOT 2000 run, and this 
spectral region is not covered in the MUSICOS 2000 run. 
In all the spectra strong emission in the Ca ii H & K 
lines and a clear emission in the He line coming from the 
primary component is observed (see Fig. U). The Ca ii 
H & K emission lines from the secondary component are 
also detected with relative wavelength shifts with respect 
to the primary in agreement with the wavelength shifts 
calculated in the Ha and the other Balmer lines. 

In our spectra the Ca ii H & K line spectral region is 
located at the end of the echellogram, where the efficiency 
of the spectrograph and the CCD decrease very rapidly, 
and therefore the S/N ratio obtained is very low, and the 
normalization of the spectra is very difficult. For these 
reasons we have not applied the spectral subtraction in 
this spectral region, and we have plotted in Fig. ^ only the 
observed spectra. The EWs measured in these spectra are 
in agreement with the strong Ca ii K emission (EW= 2.7 
A) reported by Mason et al. (1995). 

Variations of the Ca II H & K emissions with the or- 
bital phase and from one epoch to another are observed 
for the primary component, and as in the case of the Ha 
line the level of chromospheric activity is higher in the 
first observing run than in the last one. 

4.4. The Ca ii IRT lines (X8498, X8542, X8662) 

The three lines of the Ca ii infrared triplet (IRT) are in- 
cluded in our echelle spectra, except A8662 which is not in- 
cluded in the NOT 2000 observing run. In all the observed 
spectra of BK Psc a clear emission reversal is observed in 
the core of the Ca ii IRT absorption lines (see Fig. || left 
panel). After applying the spectral subtraction technique, 
only the excess emission arising from the primary compo- 
nent is detected (see Fig. || right panel). Unlike the other 
chromospheric activity indicators, no evidence of excess 
emission is detected from the cooler secondary component, 
even when in this red spectral region its relative contri- 
bution is slightly larger than in the blue one. In addition, 
the orbital phase and seasonal variations detected in the 
Ca II IRT excess emission are more smooth than in the 
other activity indicators. 

We have calculated the ratio of excess emission EW, 
■, which is also an indicator of the type of chromo- 
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Fig. 4. Observed spectra in the region of the Ca ii H & K 
and He lines. The position of these lines for the primary 
(P) and secondary (S) components are marked with short 
vertical lines 
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Fig. 5. Observed and subtracted spectra, as in Fig. ||, in the region of the Ca ii IRT (8498, 8542 A) Unes 
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solar plages, values of « 1.5-3 are measured, while 

in solar prominences the values are ~ 9, the limit of an 
optically thin emitting plasma (Chester 1991). The small 
values of the ^fsaia ratio we have found for the primary 
component in all our spectra (ranging from 1.3 to 2.4) in- 
dicate that the Ca ii IRT emission of this star arises from 
plage-like regions, in contrast to the Balmer lines that 
seem to come from prominences. This markedly different 
behaviour of the Ca ii IRT emission has also been found 
in other chromospherically active binaries (see Paper III 
and references therein). 

5. Conclusions 

In this paper a detailed spectroscopic analysis of the re- 
cently, X-ray/EUV selected chromospherically active bi- 
nary BK Psc is presented. We have analysed high resolu- 
tion cchellc spectra taken during four observing runs from 
1999 to 2001. These spectra include all the optical chro- 
mospheric activity indicators from the Ca II H & K to 
Ca II IRT lines, as well as the Li I A6707.8 line and other 
photospheric lines of interest. 

The precise radial velocities of the primary compo- 
nent of BK Psc that we have determined, by using the 
cross-correlation technique, allow us to confirm that it is 
a single-lined spectroscopic binary (SBl). In addition, as 
the chromospheric emission lines from the secondary com- 
ponent are also detected in our spectra it has been possible 
to measure the radial velocities of the secondary. We have 
used all these radial velocities to obtain, for the first time, 
the orbital solution of the system, as in the case of a SB2 
system. We have obtained a circular orbit with an orbital 
period of 2.1663 days, very close to its photometric pe- 
riod of 2.24 days (indicating synchronous rotation). By 
comparison with the spectra of reference stars of differ- 
ent spectral types we have found that the spectra of BK 
Psc is well fitted with a K5V primary component without 
any contribution of a secondary component. The mini- 
mum masses [Msirr'i) resulting from the orbital solution 
are compatible with the observed K5V primary and an 
unseen M3V secondary. As both components are main se- 
quence stars, we can classify this chromospherically active 
binary as a BY Dra system (Fekel et al. 1986). 

By using the information provided by the width of 
the cross-correlation function we have determined a pro- 
jected rotational velocity, wsini, of 17.1 km s~^ for the 
primary component. Within the errors, the radius for a 
K5V {R = 0.72Rq) is compatible with the radius obtained 
from the HiPPARCOS data {Rnip = 0.60±0.04i?o) and the 
minimum radius (i?sinz = 0.76 ± O.OSi?©) obtained from 
V sin i and the photometric period. 

The kinematics {U, V, W space-velocity components) 
and the non- detection of the Li I line in this active star 
indicate that it is an old star. 

Finally, we have analysed, using the spectral subtrac- 
tion technique, all the optical chromospheric activity indi- 
cators from the Ca ii H & K to Ca ii IRT lines. Both com- 
ponents of the binary system show high levels of chromo- 



spheric activity. In the observed spectra we have detected 
strong emission in the Ha, Ca II H & K, He, and Ca ii IRT 
lines coming from the primary component, and low emis- 
sion in the Ha and Ca ii H & K lines from the secondary. 
After applying the spectral subtraction the chromospheric 
excess emission from both components arc clearly detected 
in all the activity indicators except the Ca ii IRT lines. 
We have used a two-Gaussian fit to deblend the emission 
coming from both components. The Ha emission above 
the continuum from both components seems to be a per- 
sistent feature of this system during the period 1999 to 
2001 of our observations as well as in previous observations 
reported by other authors. In addition, the Ha emission 
of the primary exhibits a central self-absorption. The ra- 
tio that we have found indicates that the emission 
of these lines could arise from prominence-like material, 
whereas the ratio .|sm2. indicates that the Ca ii IRT emis- 
sion arises from plage-like regions. The excess emission 
EW of the Balmer lines and the Ca il H & K lines corre- 
late well with each other, but the Ca ii IRT lines show a 
different behaviour. 
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Table 6. EW of the different chromospheric activity indicators of BK Psc 
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0.46 
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0.51 



* Orbital phase calculated with the orbital period, Porb and date of conjunction, Tconj determined in this paper. 

* Data not measured due the very low S/N. 
+ Data measured with low S/N. 

^ Data for the primary and secondary components not deblended due to the orbital phase of the observation. 
^ For these lines we can observe the two components but they could not be deblended. 
^ These values have been measured at the observed spectra. 

* The He line of the primary component is blended with the Ca ll H line of the secondary. 



Table 7. Absolute surface flux of the different chromospheric activity indicators of BK Psc 

logFs in the subtracted spectrum 
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6.49/6.20 


6.19 


6.29 




NOT 00 


0.88 


* 






* 


5.41+ 


5.85/6.07 


6.44/6.48 


6.14 


6.37 




2.2 m 01 


0.17 


6.22/5.41 


6.31/" 


5.95" 


5.68/5.41 


5.72/5.02 


5.76/4.62 


6.39/6.50 


6.13 


6.23 


6.10 


2.2 m 01 


0.58 


6.50/5.34+ 


6.44/4.99+ 


5.91 


5.41/4.34 


5.58/4.98 


5.82/5.05 


6.39/6.56 


6.08 


6.22 


6.13 



Notes as in Table |^. 
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